Statistical study of short quiescent times between flares 
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Abstract. The study of the statistical distribution of short quiescent times (< 30 
minutes) between solar flares has been investigated on a ID MHD model. A power 
law behaviour is found which indicates the existence of sympathetic flaring. This 
prediction is supported by recent observations. 

Solar flares are the most violent phenomena observed in the solar 
system. It is believed that they play a fundamental role in maintain- 
ing a "hot" corona. Statistical analyses of observational data (Dennis, 
Q ■ 1985; Crosby et al., 1993; Pearce et al., 1993), in terms of probability 

distribution functions (pdfs), revealed the scale- invariance property of 
■ solar flares: the pdfs N(x) of different quantities x, such as the peak 

luminosity (P), the total energy (E) or the duration (D) of flare events, 
follow widely extended power laws, viz. N(x) ~ x~ ax (x = P, E or D) 
with the indices ap « 1.7, rs 1.5 and o>d ~ 2.2. This apparent 
£NJ \ absence of a characteristic length scale allowed us to see solar flares 

in a new way, like a Self-Organised Critical (SOC) system (Lu and 
Hamilton, 1991; Vlahos et al, 1995). Although SOC models do not 
\ make a direct connection with physical ingredients, like the magnetic 

field in the case of flares, with simple rules they are able to reproduce 
well the statistical properties observed. The limitation of SOC models 
Q-f in the predictability of flares statistics was however recently discussed 

by Boffetta et al. (1999). The authors have studied the pdf of long 
quiescent times (also called waiting times), r (for r > 6 hours), defined 
as the time intervals between two successive bursts. This observational 
analysis, based on 20 years of data obtained from the GOES sensors 
. ^ | shows clearly a power law behaviour with an index a T ~ 2.4. This kind 

^ ■ of behaviour is an indication of the existence of sympathetic flaring 

(i.e. the triggering of one flare by another) on long time-scales. It 
corroborates what we already knew for short quiescent times (Pearce 
et al., 1993; Wheatland et al, 1998). Then Boffetta et al. (1999) con- 
cluded that this result was in contradiction with SOC models of solar 
flares, where the system is driven at a constant rate (Lu and Hamilton, 
1991), and thus for which events occur randomly in time as a Poisson 
process. Note that there are many SOC models in the literature like 
nonconservative models (Christensen and Olami, 1992) which are able 
to produce power laws for the distribution of quiescent times. Therefore 
the distinction discussed above between SOC and turbulence is not so 
obvious (although it is for many other aspects!). 
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Do we really have observed a signature of sympathetic flaring at any 
time-scale ? According to Wheatland (2000) the answer is no. Indeed in 
his paper it is shown that the apparent power law behaviour observed 
for long quiescent times is simply the consequence of the wide varia- 
tion in the flaring rate during several solar cycles. However for short 
quiescent times we do have sympathetic flaring. The overabundance of 
short quiescent times (10 s to 10 minutes) constructed by Wheatland 
et al. (1998) from bursts observed in hard X-ray may be interpreted 
as a time correlation between secondary events (bursts) generated by a 
major event (flare). Then solar flares, defined as events having several 
episodes of hard X-ray emission, would be independent events. 

A first investigation of the quiescent times statistics in the context 
of magnetohydro dynamics (MHD) is described by Boffetta et al. (1999) 
through a shell model. The idea of a shell model is to reproduce the non- 
linear dynamics by including only a set of nearest-neighbour non-linear 
interactions while preserving the quadratic invariants, but otherwise 
ignoring the details of the spatial structures. In spite of its simplicity, 
this model appears to be rich enough to reproduce the observed statis- 
tics for P, D and E (see above). However, and in apparent contradiction 
with the observations, a power law behaviour is also observed for the 
pdf of the long quiescent times. 

The aim of this paper is to investigate the statistical distribution of 
short quiescent times (< 30 minutes) on a ID MHD model of coronal 
loops (Galtier and Pouquet, 1998). Several turbulent MHD models of 
coronal loops have been recently developed but few statistical studies of 
dissipative events have been reported in the literature mainly because of 
the difficulty, from a numerical point of view, to produce statistical data 
of good quality (Einaudi et al., 1996; Dmitruk et al., 1998; Georgoulis 
et al., 1998). For that reason, a reduction to a ID MHD model has been 
proposed where the magnetic loop at the origin of the flare is reduced to 
a magnetic line. Then the events generated are shocks whose dissipation 
is a source of heating for the corona. It was shown that this model was 
able to reproduce solar flares properties like the histogram of peak 
luminosity (Galtier and Pouquet, 1998; Galtier, 1999). The details and 
the results of the numerical computation considered here, with a spatial 
resolution of 4,096 grid points, are given in Galtier (1999). The main 
signal analysed in the present study is the Joule dissipation (~ (J 2 (x)), 
where J is the current density) averaged over the length of the loop for 
a simulation of 42 hours of solar activity. Each dissipative event of 
the signal defines a flare and the quiescent times are then taken to be 
the time intervals between two successive events, i.e. two successive 
maxima, following the criterion also used in Wheatland et al. (1998) 
and in Boffetta et al. (1999). The computation of the histogram of 
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Figure 1. Histogram JV(r) of the short quiescent times (in second) between two 
successive dissipative events. The straight line corresponds to a slope of —2.3. 

the quiescent times is displayed in Figure 1 in log-log coordinates; the 
unit time chosen is the second. A clear power law behaviour is revealed 
over a range of time from 3 minutes to 30 minutes, with an index 
a T f» 2.3. It indicates the existence of sympathetic flaring in agreement 
with the observations made by Pearce et al. (1993) for a range of short 
quiescent times from 1 to 60 minutes, and with the observations made 
by Wheatland et al. (1998) where an overabundance of short quiescent 
times is seen from 10 s to 10 minutes. In fact in the latter case the data 
can also be fitted with a power from 10 2 s to 10 4 s. We have to emphasise 
here that, (a) the power law index obtained with the MHD simulation 
is different from the observations. However it is claimed that the most 
important result is the power law behaviour which implies the existence 
of sympathetic flaring for short time-scales, (b) Some precautions have 
to be taken with the result obtained by Pearce et al. (1993) since the 
authors did not take into account the variation in the flaring rate like 
in Wheatland et al. (1998). 

The physical picture for long and short quiescent times appears to 
be very different. For long quiescent times (longer than the lifetime of 
a flare) no time correlation seems to be observed, however for short 
quiescent times (of the order or smaller than the lifetime of a flare) we 
do observe a time correlation. As it is proposed by Wheatland et al. 
(1998) the origin of sympathetic flaring in the MHD model is possibly 
due to the correlation between secondary events (bursts) generated 
by a major event (flare). In fact typical dissipative events last several 
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minutes (Galtier, 1999) and are often composed in a main flare followed 
by secondary events. The detailed analysis, including spatial informa- 
tion along the coronal loop, on typical examples of dissipative events 
reveals indeed a non-trivial correlation between bursts (Galtier, 2000) 
for short time-scales. This spatial information available in the MHD 
model appears to be an important ingredient to understand the origin 
of sympathetic flaring. 
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